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Theory  of  linear  optical  properties  of  fractal  clusters  is  developed^  The  theory  is 
based  upon  the  exact  properties  of  dipole  polarizability  and  assumption  of  the  existence 
of  scaling  for  the  dipole  excitations  (eigenstates)  of  the  fractal.  This  assumption  is  self- 
consistentlv  validated  by  the  results  of  the  theory  and  is  also  confirmed  by  numerical 
-timulation  in  the  framework  of  the  Monte-Carlo  method.  Using  exact  relations  and  the 
scaling  requirements,  it  is  shown  that  the  fractal  absorption  and  density  of  eigenstates  scale 
with  the  same  exponent  d0  —  1.  .The  index  d0,  which  is  called  the  optical  spectral  dimension, 
is  the  counterpart  of  the  Alexander-Orbach  fracton  dimension  d.  The  dispersion  law  for 
the  fractal  eigenstates  is  found:  it  is  governed  by  d„  and  the  Hausdorff  dimension  D.  Using 
this  law.  the  condition  of  the  scaling  is  determined.  According  to  this,  the  scaling  region 
occupies  the  center  of  the  absorption  band  of  the  fractal.  The  spectral  dimension  d0  is 
found  for  three  types  of  the  fractals  studied  from  the  results  of  the  numerical  simulation. 
The  wings  of  the  fractal  absorption  contour  are  described  by  introducing  the  model  of 
diluted  fractals  with  the  use  of  the  binary  approximation.  The  general  conclusion  is  that 
large  fluctuations  of  local  fields  dominate  in  determining  the  fractal  optical  spectra. 

t  Also  with  the  Institute  of  Automation  and  Electrometry,  Siberian  Branch  of  the 
USSR  Academy  of  Sciences.  630090  Novosibirsk.  USSR 


1.  Introduction 


Fractals,  as  introduced  by  Benoit  Mandelbrot  over  ten  years  ago.  are  scale  self- 
similar  mathematical  objects  possessing  nontrivial  geometrical  properties,il.2!.  TThere  exist 
various  physical  realizations  of  fractals  J3 .41.  and  here  we  shall  consider  what  we  believe 
to  be  one  of  the  most  important  such  realizations,  namely,  fractal  clusters.  Attention  will 
be  paid  mainly  to  their  optical  properties. 

"  A  fractal  cluster  is  a  system  of  interacting  material  particles  called  monomers. 
This  system  is  self-similar  (in  a  statistical  sensei  with  respect  to  scale  transformations  in 
an  intermediate  region  of  sizes  r.  /?,,  <<  ;•  Rr.  where  R  is  a  characteristic  separation 
between  nearest  monomers,  and  /?„  is  the  cluster  total  size.  For  the  sake  of  brevity,  such 
(dusters  will  simply  lie  called  fractals,  since  other  realizations  of  fractals  will  not  be  used 
in  this  paper. 

Fractals  are  widespread  in  nature  [4.5].  Products  of  a  broad  class  of  diffusion- 
controlled  aggregation  reactions  in  solutions  and  gases  can  be  labelled  as  fractals.  Accord¬ 
ingly.  fractals  are  particles  in  colloidal  solutions,  sols  and  gels,  and  soot  and  smoke,  and 
most  macromolecules  are  fractals.  Rough  surfaces,  disordered  layers  on  surfaces,  porous 
objects  ( microfilters  including  activated  coal,  porous  glasses,  aerogels,  many  heterogeneous 
catalysts,  etc.  i  also  possess  fractal  structure  in  an  intermediate  ranee  of  sizes. 

A  consequence  of  the  scale  self-similarity  of  fractals  is  a  power  dependence  of 
correlation  functions  on  coordinates.  In  particular,  the  pair  (density-density)  correlation 
function  in  the  intermediate  region  has  the  scaling  form 

D  (  r  \D-3 

,,r|  =  IojjU)  (1) 

where  the  index  D  is  called  the  fractal  (external  Haussdorff)  dimension.  Equation  (1)  is 
also  a  definition  of  the  constant  Rq .  A  consequence  of  ( 1 )  is  the  scaling  behavior  of  the 
number  .V  of  monomers  in  a  fractal  and  its  mean  density  p: 

-V~t  RJRn)°  .  «Rc/f?0)D-3  .  (2) 

A  fractal  is  called  nontrivial  if  D  is  less  than  the  dimension  of  the  space  it  is  embedded 
in.  or  in  our  case  of  a  cluster,  if  D  <  3.  As  one  can  see  from  (2),  the  mean  density  of  a 


nontrivial  fiactal  is  asymptotically  i  for  Rr  Rn)  zero.  This  feature,  together  with  strong 
pair  correlation  (1).  is  the  reason  for  the  great  role  that  tiuctuations  play  in  a  nontrivial 
fractal. 

Among  the  physical  properties  of  fractals,  the  optical  properties  have  been  the 
least  studied.  For  example,  in  a  set  of  proceedings  [4]  several  years  ago  devoted  to  the 
physical  properties  of  fractals,  there  are  no  papers  on  the  optical  properties.  The  linear 
optical  (dipole)  polarizabilities  of  fractals  have  been  studied  theoretically  [6-9].  However, 
in  Refs.  16. 7j  based  on  various  modifications  of  the  mean-field  method,  strong  fluctuations 
i  of  density,  local  fields,  etc.)  in  fractals  are  not  completely  taken  into  account.  In  later 
papers  S.9].  fluctuations  in  the  immediate  surroundings  of  a  monomer  are  seen  to  play  a 
decisive  role.  There,  the  binary  approximation  was  employed:  the  interaction  of  a  monomer 
with  only  its  nearest  neighbor  was  accurately  taken  into  account,  while  the  effect  of  other 
monomers  was  simulated  by  the  Lorentz  field. 

If  monomers  are  high-quality  optical  resonators,  for  example,  particles  of  coinage 
ie.g..  Ag  and  Au)  metals  or  other  metals  possessing  strong  luster,  strongly-fluctuating  lo¬ 
cal  fields  in  a  fractal  can  significantly  exceed  the  exciting  external  field.  This  very  feature 
brings  about  a  giant  enhancement  of  Raman  scattering  from  fractals  [S.9].  Nonlinearities, 
as  usual,  increase  the  effect  of  fluctuations  and  related  to  this  are  the  huge  magnitudes 
of  noniinear  polarizabilities  of  fractals  predicted  in  [10.11].  This  prediction  was  then  con¬ 
firmed  experimentally  [12]  for  optical  phase  conjugation  in  a  degenerate  four- wave  mixing 
process  in  fractal  silver  clusters.  Photomodification  of  such  clusters,  which  is  selective  in 
i  he  wavelength  and  polarization  of  the  radiation,  was  later  observed  [13].  These  estab¬ 
lished  effects  are  rich  in  physics  and  promise  interesting  applications.  Therefore,  they  are 
worthwhile  for  more  detailed  theoretical  study. 

The  aim  of  this  paper  is  to  develop  a  theory  of  optical  properties  of  fractals,  taking 
full  account  of  their  fluctuation  nature  without  the  restriction  of  the  binary  approximation. 
The  theory  is  based  on  a  scale-invariant  approach.  In  Sec.  2  we  present  the  basic  equations, 
and  in  Sec.  3  we  examine  the  polarizability  in  the  collective  region.  The  fractal  optical 
response  at  the  wings  of  the  absorption  contour  is  addressed  in  Sec.  4.  results  of  numerical 
simulation  are  given  in  Sec.  5.  and  Sec.  6  provides  a  concluding  discussion. 
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2.  Basic  equations  and  exact  properties  of  solutions 


Let  us  consider  a  model  of  a  cluster  as  a  system  consisting  of  .V  polarizable  particle 
called  monomers)  located  at  the  points  r,.  i  =  The  total  size  of  the  cluster 

is  supposed  to  be  much  less  than  the  wavelength.  Therefore,  the  interaction  between 
monomers  at  optical  frequencies  is  principally  the  interaction  of  their  transitional  multipole 
moments.  The  main  part  of  this  work  deals  with  interactions  at  distances  which  are  much 
greater  than  the  monomer  size,  which  is  valid  for  both  the  collective  region  (see  Sec.  3) 
and  diluted  fractals  in  the  binary  region  (see  Sec.  4).  In  this  case  the  most  important  part 
of  the  multipole  interaction  is  the  dipole  one.  However,  if  the  monomers  are  conducting 
particles,  the  charge  transfer  between  them  i ohmic  current)  should  also  be  considered.  In 
Maxwell's  equations,  the  dipole  interaction  in  a  continuous  medium  is  described  by  the 
displacement  current  jj  =  >  e/c )(d'E/dt)  =  ie\uj/c)E.  where  e  is  the  dielectric  constant.  E 
is  electric  field.  is  the  optical  frequency.  The  ohmic  current  jn  =  crE,  where  is  a  is  the 
conductivity.  The  ratio  \ju/Jd\  shows  us  how  important  the  ohmic  coupling  is  with  respect 
to  the  dipolar  one.  Using  the  above  expressions,  in  which  e  is  taken  to  be  the  dielectric 
permittivity  of  the  monomer  material,  to  describe  the  pair  interaction  between  nearest 
monomers,  we  can  obtain  the  upper  estimate  for  \jn/jd\-  The  reason  why  this  estimate  is 
the  upper  one  is  that  the  ohmic  contact  between  monomers  may  have  high  resistance  due 
to  its  small  area,  surface  oxidation,  spatial  separation,  etc.,  while  the  dipole  interaction  is 
determined  by  geometrical  factors  and.  therefore,  is  the  same  order  of  magnitude  as  in  the 
continuous  medium.  Using  the  familiar  relation  a  —  (j,'/47rllme.  we  obtain  the  estimate 
jq l jd j  <  |Ime/Ree|.  In  fact,  in  the  visible  portion  of  spectrum  this  ratio  is  small  for  most 
metals  since  —  Ree  is  large,  and  for  semiconductors  and  dielectrics  since  line  is  small.  For 
example,  for  A  =  650  nm.  \jn/jd I  =  0.07  and  0.02  for  Ag  and  Si.  respectively  (the  optical 
data  are  taken  from  Ref.  [14]).  We  shall  consider  materials  of  the  types  indicated  above 
and.  therefore,  shall  neglect  the  ohmic  coupling. 


The  dipole  (transition!  moments  induced  on  monomers  obey  the  well-known  sys¬ 
tem  of  equations 

.v 

-  r  ...  ,:_>1 

ij  •  '3) 


=  \0 E™  -  xo  Y'  krf  -  3n^)n(;j)]  r~\iu 


j~  i 


■J 


where  the  Greek  subscripts  stand  for  tensor  components  i  summation  over  repeated  indices 
is  implied),  and  the  Latin  indices  stand  for  ordinal  numbers  of  monomers:  \o  is  the  dipolar 
polarizability  of  an  individual  ( isolated)  monomer,  taken  to  be  isotropic  to  avoid  inessential 
complications;  Et0)  is  the  electric  rieid  of  the  external  (exciting;  radiation  at  the  site  of 
the  ith  monomer:  ri;  =  r,  —  r,.  and  n^’  =  r 

In  the  electrodynamics  of  continuous  media.  Eq.  (3)  is  treated  in  the  mean-field 
approximation,  in  which  it  is  reduced  to  the  integral  equation  for  the  mean  dipole  mo¬ 
mentum  d.  d  =  \oE|0)  -  \o  /  :d  -  3r(rd)/r-lr-^</;r)dr.  Taking  (1)  into  account,  one 
can  easily  see  that  for  the  nontrivial  fractal  the  right-hand  side  of  this  equation  diverges 
at  smail  distances  as  rD~'.  It  means  that  the  most  important  is  the  contribution  of  the 
nearest  monomers.  This  feature  was  exploited  in  recent  work  ;$-llj  in  which  the  binary 
approximation  has  been  developed  to  exactly  take  into  account  the  field  of  the  closest 
monomer,  with  the  rest  of  monomers  treated  by  the  mean  field.  Due  to  the  divergence  at 
r  — *  0.  this  approach  is  sensitive  to  the  small-scale  structure  of  the  fractal,  and  the  po¬ 
larizability  found  does  not  obey  scaling  in  the  general  case.  In  the  present  paper  (Sec.  4), 
we  introduce  the  model  of  diluted  fractals,  which  have  the  simplest  possible  small-scale 
structure  and  obey  scaling  both  at  intermediate  and  small  scales.  In  this  model,  the  binary 
approximation  is  shown  to  correctly  describe  the  wings  of  the  absorption  contour  of  the 
fractals. 

We  show  below  i  Sec.  3)  that  the  central  part  of  the  spectra  can  self-consistently  be 
described  as  scale-invariant.  This  means  that  the  polarizability  is  determined  by  collective 
i  long-scale  correlated)  fluctuations  of  the  local  dipole  momenta  d,.  Since  these  fluctuations 
dominate,  the  the  mean-field  approach  is  not  applicable.  In  this  respect,  the  fractals 
resemble  systems  near  the  phase  transition  point.  In  the  collective  region,  the  fractal 
optical  properties  are  not  sensitive  to  the  small-scale  structure,  scaling  takes  place,  and 
the  diluted  fractals  has  the  same  polarizabilities  as  the  original  ones. 

To  develop  the  general  theory  of  the  optical  properties  of  fractals,  we  shall  use  the 
approach  based  on  the  expansion  of  polarizability  in  terms  of  the  exact  eigenvectors  and 
eigenstates  of  the  interaction  operator  (this  is  similar  to  the  Lehman  expansion  in  field 
theory;.  Also,  such  an  expansion  turns  out  to  be  the  most  efficient  computing  method  for 
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carrying  out  our  numerical  simulations  in  Sec.  5. 

To  perform  the  above  indicated  expansion,  let  us  represent  (3)  in  the  form  of  an 
operator  equation  by  introducing  the  matrix  U'  with  elements 

(laltrl jj)  =  (  i  *  j:  (4) 

l  0.  i  =  J  • 

which  operates  in  a  3Ar-dimensional  vector  space  of  the  vectors  d  and  E  with  components 


( mid)  =  dlQ  .  (miff)  =  Exc 


(5) 


Then  the  basic  system  of  equations  (3)  acquires  the  form 


Z  +  l V)d  =  E[0)  . 


(6) 


where  the  complex  variable  Z  with  real  (  —  X)  and  imaginary  (—<!>)  parts  is  introduced. 

Z  =  -(X  +  i6)  =  xo1  ■  (?) 

The  symmetric  real  matrix  W  is  diagonalized  by  the  orthogonal  transformation 

i  UTt  =diag(u?n)  .  UUT  =  1  .  (S) 

where  diag(ien),  is  the  diagonalized  matrix  consisting  of  eigenvalues  wn,  and  the  super¬ 
script  T  denotes  the  transpose.  The  formal  solution  of  the  basic  equation  (6)  has  the 
form 

d=rTdiag(— - - )UEi0)  .  (9) 

Z  4-  wn 

Expressing  the  required  quantities  d,Q  in  terms  of  the  eigenvalues  wn  and  components  of 
the  corresponding  eigenvectors  (n|m).  one  obtains  from  (9) 

-/,a  =  ^(n|T'!ta)(n|tf|jJ)(Z  +  u.n)-1£,^)  .  (10) 


We  consider  below  the  clusters  whose  size  Rr  is  much  smaller  than  the  radiation 
wavelength  A.  Then  the  external  field  E^J  may  be  considered  homogeneous,  i.e..  not 


iepending  on  j.  In  this  case  the  liipole  moment  induced  on  the  /lh  monomer  is  expressed 


in  terms  of  the  corresponding  linear  polarizability  \^j  as 


I  t  )  U  ) 


It  is  clear  from  (101  that  \'Jd  is  given  by 


ill) 


\'^j  =  (n\U\ia)(n\U\j  3)  {Z  -r  wn)  (12) 

The  polarizability  as  an  analytical  function  of  Z  differs  from  that  expressed  in  terms  of  fre¬ 
quency  In  particular,  this  function  is  not  "real"  in  the  sense  that  V.'jV-Y)  =  \ ^ j (  —  -Y  'i . 
Nevertheless,  it  can  be  shown  directly  from  the  solution  i  12)  that  the  conventional 
Kramers- Kromg  formula  in  A'  is  valid. 


Re\ 


1 1) 
rt  J 


(13) 


where  P  denotes  the  principal  value  of  the  integral.  We  shall  obtain  also  some  other  exact 
relations  for  \^j  which  are  a  consequence  of  the  fact  that  the  eigenvalues  wn  are  real,  the 
matrix  U  is  orthogonal,  and  6  is  positive  corresponding  to  normal  absorption  (we  assume 
no  population  inversion  i. 


From  the  form  of  the  solution  i  12).  the  exact  sum  rules  follow: 

X 

3  j  hn\['t)j(X)dX  =  .  P  J  Rey^A'idA'  =  0  .  (14) 

—  TO  —  -^C 

From  I  14)  it  can  be  seen  that  the  absorption  integral  is  conserved:  it  is  the  same  as  for  an 
isolated  monomer.  The  absorption  cross  section  differs  from  Im\ ad  only  by  a  trivial  factor 
4 ~k.  where  k  is  the  wave  vector.  As  a  characteristic  of  absorption  we  shall  use  Im\a'd- 

We  now  obtain  an  exact  relation  which  is  a  counterpart  of  the  optical  theorem. 
For  this  purpose,  from  (  12)  and  taking  (3)  into  account,  we  find  a  quadratic  form  of  the 
polarizabilities: 


For  definiteness,  let  us  suppose  that  the  external  field  E  is  directed  along  the  r-axis.  From 
IS)  and  1 19)  we  obtain 

JjE1|7lE",)j2  =d(l-t-A'-/d‘2)Im\„  .  (20) 

I 

This  relation  resembles  the  Callen- Welton  formula.  Earlier  a  result  had  been  obtained 
with  the  use  of  the  binary  approximation  (see  Eq.  (45)  in  Ref.  (9])  which  differs  from  (20) 
only  by  some  notations.  This  result  was  the  basis  to  develop  a  theory  of  enhanced  Raman 
scattering  from  fractals. 

We  emphasize  that  the  above  relations  ( 13).  ( 14).  ( IS)  and  (  20)  are  exact.  They  are 
valid  for  an  arbitrary  cluster  ( fractal  or  not)  without  any  averaging  over  cluster  ensemble  or 
orientation.  It  is  only  essential  that  the  cluster  consists  of  monomers  with  pair  interactions 
between  them  which  give  a  linear  response  to  an  external  field  and  possess  normal  (not 
population-inverted)  absorption. 
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Let  us  now  invoke  rotational  symmetry.  After  averaging  over  the  orientation  of 
a  cluster  as  a  whoie  i  denoted  beiow  as  the  polarizability  tensor  is  reduced  to  a 


'cater. 


\aJ/  —  \Taii  •  \  —  ^  VorQ 

For  the  first  moment  of  \.  we  ran  find  directly  from  (  12) 


XI  ?C 

P  j  ATm\i  A'hLY  =  P  J  Rex(X)dX -t 


\ia)wn(n\U'  j  a). 


Accordingly  to  i  14  u  the  first  t  rn  on  the  right  side  of  (22)  equals  zero.  The  second 
‘erm  with  the  iieip  of  (S'  is  reduced  to  the  form  [in  IF  ja  u  Using  the  fact  that  the 
interaction  tensor  (41  becomes  zero  when  convoluted  over  the  indices  o  and  .i.  one  comes 
to  the  conclusion  that  the  second  term  is  also  equal  to  zero.  In  such  a  way.  we  obtain  the 


exact  sum  rule 


P  /  ATm\'.  X  FIX  =0  . 


meaning  that  the  absorption  center  of  a  cluster  is  at  X  =  0. 

The  universality  of  the  above  obtained  relations  is.  in  particular,  due  to  the  use  of 
natural  variables  A"  and  <"  for  the  present  problem.  These  variables,  indeed,  depend  upon 
Tie  frequency  w.  but  nor  in  a  universal  way.  In  what  follows,  it  is  important  to  establish 
the  region  in  which  the  spectral  variable  A  is  changing.  Let  us  consider  the  dependence 
X ( w  i  in  two  different  models. 

First,  let  a  monomer  possess  a  single  isolated  resonance  with  frequency  •  relax¬ 
ation  rate  (homogeneous  width,  f  and  transition  dipole  matrix  element  d  12-  In  this  case 


where  Q  -  —  w() .  One  can  see  that  A  is  a  relative  detuning  and  b  is  a  relative  width  of  the 

resonance  in  an  isolated  monomer.  We  point  out  that  the  first  of  the  relations  (  14  )  together 
with  (23)  reproduces  the  well-known  dipole  sum  rule  for  the  polarizability.  Note  that  the 
minimum  value  of  T  is  equal  to  the  natural  (radiative)  linewidth  Tn  =  4fc3 jdi2j"/3fi. 
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Substituting  F  =  1  „  into  Eq.  i  24).  we  obtain  tor  the  resonance  i  *2  =  ill  riie  upper  estimate 
I:n\D  \!.  which  is  much  greater  than  the  volume  of  the  monomer,  in  ai^ord  with  the 
quantum-mechanicai  Breit-W  igner  iormuia.  In  real  clusters,  rhe  width  1’  is  determined 
uy  radiativeless  dissipation  and  dephasing,  and  T  ">  F .  Thus,  the  magnitude  of  Im\o  is 
mucii  smaller  than  A5,  though  it  may  still  be  much  greater  than  the  monomer  size  due  to 
resonant  enhancement  i  cf.  ‘he  example  below). 

The  second  modei.  which  is  realistic,  in  particular,  for  clusters  obtained  bv  aggre¬ 
gation  m  colloidal  solutions,  considers  rhe  monomer  as  a  macroscopic  sphere  with  radius 
R .  consisting  of  rhe  material  with  a  complex  dielectric  permittivity  ;  —  s'  —  is".  The 
•  >i uarizability  ot  su-di  a  sphere  is  given  by  rhe  weil-known  expression 

V„  =  R3m(s  -  l)(c  4-2)-'  .  (25) 

For  the  metallic  sphere,  the  permittivity  is  well  described  by  the  Drude  formula 

£  =  '-o  ~  ^p/[w.’(u,’  4-  i~f  j]  .  (26) 

where  fo— const  is  rhe  contribution  of  the  interband  transitions  (e.g..  for  silver  co  ~  5).  ujv 
rilt’  electron  plasma  frequency.  and  -  is  the  electron  collision  frequency.  For  most  metals. 
:u  particular.  lor  the  noble  metals.  *.  <$2  and.  consequently,  s"  ■'<  .  In  this  case,  from 

1  1  )  !‘U<1  (25).  one  can  see  that  X l )  turns  to  zero  at  the  point  of  the  surface  plasmon 
resonance.  _•  =  j,.  where  s'  —  —  2.  From  (26)  it  follows  that  w,  mp(co  -r  2)-1^*.  In  the 
vicinity  of  the  surface  plasmon  resonance,  the  dependence  A’(dci  reduces  to  (24)  with  the 
following  parameter  values:  =  w  -  w,.  (d12 j‘  =  3R3mh^p/  [2(c0  -r  2)J/Jl.  T  =  7/2.  Also, 

it  follows  from  (25)  and  (2C)  that  for  w  <c  the  variable  X  tends  to  -1.  In  the  vicin 
of  the  frequency  the  variable  Z  (7)  experiences  strong  dispersion 


Z  ~  jRjiso  -  l)  / [- 


-+-  i-v/21  . 


It  follows  from  1  27)  that  in  a  small  vicinity  of  u the  variable  A  changes  sign,  acquiring  a 
large  absolute  value  iimited  by  A'mar  =  7- 
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3.  Polarizability  of  a  fractal  in  the  collective  region 

To  obtain  the  above  exact  results  we  did  not  use  any  specificity  of  fractals.  Those 
resuits,  however,  will  now  be  employed  to  obtain  a  description  of  the  collective  polarizabil¬ 
ity  of  fractals.  Self-similarity,  which  is  a  fundamental  property  of  fractality.  means  that 
a  fractal  reproduces  itself  when  the  spatial  scale  is  changed.  Since  in  general  a  fractal 
is  a  random  object,  this  reproduction  has  statistical  meaning.  A  change  of  length  scale 
brings  about  a  simultaneous  change  of  the  scale  of  eigenvalues  wn  and  the  variable  X .  As 
a  consequence  of  this,  intermediate  asymptotic  values  of  observed  quantities  as  functions 
of  ir  and  X  should  be  power-like  (scaling). 

We  shall  describe  first  Im\(A"  i  1  absorption  i  and  later  consider  Re\  > A* )  with  the 
aid  of  *  13).  From  (12)  and  taking  i7)  into  account,  we  see  that  the  universal  scaling 
behavior  of  Im\(A’)  may  exist  only  if  A’  A>  6.  Then 

Im\i  X)  =  —  ^ / (n|£/|ia)  (n\U\ja)6(X  -  wn)^  .  (28) 

ijn 

where  •••/  denotes  an  average  over  the  ensemble  of  fractals.  From  (28)  it  follows  that 
absorption  is  determined  by  those  eigenstates  for  which  wn  —  X . 

The  eigenstates  of  any  large  system,  including  the  fractal,  are  characterized  by  a 
coherence  length  (which  turns  out  to  be  the  wavelength  in  the  trivial  case  D  =  3).  within 
which  the  excitations  of  monomers  are  strongly  correlated.  This  coherence  length  is  a 
function  of  the  generalized  frequency  A’,  which  we  shall  denote  as  L  \ .  The  eigenstates  are 
collective,  i.e..  delocalized  over  many  monomers,  if  L  \  7>  f?r>.  The  scaling  should  exist  in 
the  intermediate  region 

Rn  «  Lx  <  Rr  ■  (29) 

Our  basic  assumption  is  that  the  dependence  of  L\  on  A’  allows  the  existence  of  a  finite 
interval  of  X  in  which  (29)  is  satisfied.  Assuming  so.  we  shall  develop  the  scale-invariant 
theory  of  the  collective  excitations  of  fractals.  This  theory  will  eventually  allow  us  to  find 
the  dependence  of  L  \  from  the  scaling  properties  of  relevant  quantities  and.  in  such  a  way. 
to  seif-consistentlv  confirm  the  basic  assumption.  It  is  also  supported  by  the  results  of  the 
numerical  simulation  (Sec.  5). 


In  the  sealing  region  i29).  the  form  of  the  absorption  contour  Im \ l A  )  should 
be  power-like.  The  exact  symmetry  of  the  eigenstate  problem  of  the  U  operator  with 
respect  to  a  sign  change  of  w„  is  absent.  Consequently,  parity  in  X  is  not  held,  suggesting 
that  the  indices  for  X  >  0  and  A’  <  0  could  be  different.  However,  the  sum  rule  (23) 
together  with  the  self-similarity  of  a  fractal  guarantees  that  these  indices  coincide,  and  the 
corresponding  coefficients  also  coincide.  Thus,  in  the  region  (29)  (provided  jA'j  S),  a 
symmetric  scaling  dependence  of  Imy  is  valid,  with  coefficients  which  are  determined  by 
Ro  and  can  be  estimated  from  dimensionality  arguments. 

lm\(X)  ^  R3r){R*\X\)da~l  .  (30) 

where  d0  is  an  index  which  we  call  the  optical  spectral  dimension.  From  the  convergencv 
requirement  for  the  first  of  the  sum  rules  !  14),  it  follows  that  d0  >  0.  We  should  emphasize 
that  the  symmetry  of  absorption  in  X  is  not  an  exact  property,  but  is  valid  only  in  the 
collective  region,  and  is  violated  in  the  binary  region  (see  Sec.  4) 

Let  us  consider  the  transformation  properties  with  respect  to  a  change  in  the 
minimum  size  Ro  (renorm-transformation).  We  take  the  initial  value  as  Ro  =  1.  i.e..  we 
measure  lengths  in  the  unit  of  the  initial  size  Rq  .  Increasing  Ro  means  coarsing  of  the 
spatial  precision  with  which  a  fractal  is  viewed.  Such  coarsing  can  be  achieved  in  the 
following  manner  [loj.  We  first  isolate  in  a  fractal  fluctuations  with  sizes  on  the  scale  of 
/  1.  We  then  consider  monomers,  which  form  such  a  fluctuation,  as  a  new  composite 

(renormalized)  monomer.  From  the  property  of  self-similarity,  it  follows  that  the  obtained 
cluster,  consisting  of  renormalized  monomers,  is  a  fractal  with  the  same  critical  exponents, 
and  Ro  ~  (■  The  renorm-transformation  conserves  the  form  of  the  absorption  contour  (30), 
provided  the  collectivity  criterion  (29)  is  met. 

We  now  determine  the  renormalization  law  for  A’  and  find  the  coherence  length 
L.\  ■  The  renorm- transformation  should  conserve  the  total  absorption  of  a  fractal.  This 
requirement,  taken  together  with  (2)  and  (30).  reduces  to 
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where  cue  proportionality  is  understood  in  the  sense  of  the  dependence  on  f?0  land  not  on 
Rr ).  From  (31)  it  is  clear  that  A'  transforms  in  the  following  way: 

Y'  ^  — t3d„  — D)/(d„  —  1  )  (  30) 

We  note  that  in  the  case  of  a  trivial  fractal  ( D  =  3).  j  A’  j  oc  R ^ !  independently  of  the  index 
dr, ,  i.e..  we  have  the  well-known  result  Re\o  a  R^,  which  demonstrates  the  proportionality 
of  the  polarizability  of  a  composite  monomer  to  its  volume.  We  should  point  out  that  the 
transformation  law  (32)  is  derived  from  general  principles,  and  the  single  condition  of  its 
validity  is  the  existence  of  the  scaling  behavior  of  absorption  1  see  below). 

For  strongly-disordered  systems  such  as  fractals.  Alexander  lot  has  assumed  that 
collective  states  with  a  given  frequency  (whose  role  in  this  case  is  played  by  the  variable 
A"  ),  are  characterized  by  a  single  coherence  length  which  is  an  analog  of.  simultaneously, 
the  wavelength  and  localization  radius.  This  property,  called  strong  localization,  is  a 
consequence  of  the  self-similarity  of  a  fractal  and  its  collective  states.  For  a  trivial  (three- 
dimensional)  system,  this  property  in  the  general  case  can  obviously  be  violated:  the 
excitation  localization  radius  can  essentially  exceed  the  wavelength,  which  is  the  situation 
of  Anderson  localization. 

The  coherence  length  L \ .  as  other  observable  quantities  in  the  scaling  regime, 
should  net  depend  on  both  the  maximum  Rc  and  minimum  length  Ro.  The  latter  require¬ 
ment  means  that  L,\  should  be  an  invariant  of  the  renorm-transformation.  Taking  into 
account  Eq.  (32).  one  can  make  sure  that  it  possible  to  construct  only  one  quantity  which 
possesses  the  above  indicated  properties  and  has  the  dimensionality  of  length. 

Lx  x  R„(^|W|)U’'1>/(3"D1  .  (33) 

Thus,  for  the  present  problem.  Alexander’s  strong  localization  follows  from  scaling.  The 
relation  reciprocal  to  (33)  is  the  dispersion  law  for  the  excitations  of  the  fractal: 

\X(Lx)\  *  R^  (Lx/Rof~D)nd°~l)  ■  (34) 

This  dispersion  law  for  the  collective  dipole  excitations  of  the  fractal  can  be  considered  as 
a  theorem  based  upon  the  general  relations  and  the  scaling  assumption. 
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For  the  trivial  system  i  D  =  3).  the  expression  (331  is  singular,  which  is  not  sur¬ 
prising  (see  above).  Nevertheless,  this  singularity  has  simple  physical  meaning,  and  to  see 
this,  we  consider  the  dispersion  law  (34).  For  D  —  3.  it  follows  from  (34)  that  X(L\) 
does  not  depend  on  L  y.  he.,  on  the  excitation  wavelength.  jA'I  ~  ] .  This  estimate 

exactly  corresponds  to  the  dispersion  law  ( existence  of  a  spectral  gap )  of  long- wave  surface 
plasmons.  Note  that  in  the  trivial  limit  (D  — *  3)  the  dipole  excitations  of  a  fractal,  which 
give  contribution  to  the  optical  absorption,  tend  exactly  to  the  surface  plasmons  since  the 
radius  of  the  fractal  Rr  is  assumed  to  be  much  less  than  the  radiation  wavelength  A.  To 
demonstrate  this,  consider  an  averaged  trivial  cluster  as  a  microscopic  sphere,  for  which 
the  surface  plasmon  resonance  condition  is  e'  =  —2  (cf.  discussion  of  Eqs.  (  25)-(27)),  where 
e  is  the  dielectric  permittivity  of  the  sphere  material.  Recalling  that  e  =  1  4-  4-f?^  *\o 
and  taking  (7)  into  account,  we  obtain  for  the  resonance  X  =  —  ^R^* .  in  accord  with  the 
above  estimate. 

For  a  normal  dispersion  law.  the  "frequency''  X  should  decrease  with  the  wave 
vector,  i.e..  with  an  increase  of  L_\-  Taking  into  account  (34),  this  leads  to  the  limitation 
do  <  1.  which  we  accept  (as  will  be  confirmed  later  by  numerical  simulations).  We  shall 
not  use  the  renorm-transformations  below  and  shall  return  to  the  original  value  Rq  =  1. 
Then  the  condition  (29).  which  determines  the  scaling  region,  together  with  (2)  and  (33) 
acquires  the  form 

_V(3/D—  n/(i—  do)  <  x  <  x  (35) 

Hence,  this  region  is  spread  to  small  |  _V  j  and,  consequently,  occupies  the  center  of  the 
absorption  band  (cf.  (24)).  Obviously,  the  necessary  condition  for  scaling,  |A’|  6.  is 

compatible  with  (35)  only  if  the  Q-factor  of  the  resonance  in  a  monomer  is  high: 

Q  =  1/6  >  1  .  (36) 

As  follows  from  (35),  when  a  fractal  tends  to  a  trivial  one  (D  —*  3).  the  region  of  collective 
behavior  in  A  degenerates  to  a  point  corresponding  to  the  frequency  of  surface  plasmons 
(see  above).  Thus,  a  finite  size  of  the  band  of  collective  states  of  a  fractal  is  entirely  due 
to  its  nontriviality. 
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The  scaling  region  <  351  can  be  physically  realizable  only  if  the  interval  of  values  of 
the  function  A'(-’ )  (  7)  includes  it.  This  interval  depends  on  which  monomers  the  cluster  is 
built  with.  Let  us  discuss  in  this  connection  the  above  considered  model  of  the  spherical 
monomer  ( see  |  25)-(27)  and  the  corresponding  text ).  The  monomer  radius  Rm  is  supposed 
to  be  on  the  order  of  Rr>  =  1  (in  the  present  units).  There  exist  small  values  of  A”,  including 
zero,  near  the  frequency  ^  if  the  surface  plasmon  resonance  is  well  pronounced,  the 

condition  of  which  7  <C  coincides  with  (36).  In  practice,  it  is  well  met  for  many  metals, 
e.g..  for  silver  y/u,  <  1/30.  with  the  resonant  frequency  corresponding  to  the  near-UV 
wavelength  of  approximately  300  nm.  From  the  discussion  of  Eqs.  ( 25 )-( 27 ).  it  follows  that 
in  the  most  interesting  band  of  A  >  300  (i.e..  soft-UV.  visible  and  IR  light  I  the  variable 
A’  belongs  to  the  interval  from  -i  to  0.  which  includes  the  whole  scaling  region  i35).  The 
variable  A'  is  in  the  interval  from  0  to  a  large  value  of  ~  u;p/y  when  w  is  from  u,q  to  =  uj 
ifar-UV  band),  and  large  negative  values  of  A'  occur  for  >  ~'d- 

To  summarize  the  above  obtain  results,  let  us  point  out  that  we  have  initially 
assumed  there  to  exist  the  region  (29)  of  the  collective  fractal  excitations.  Starting  from 
this  point  and  using  exact  relations,  we  have  shown  that  the  absorption  is  given  by  the 
scaling  law  (30).  found  the  scaling  properties  (34)  of  the  variable  Ah  and  determined  the 
dispersion  law  for  the  fractal  excitations  (33)  (or  (34)).  From  this  law.  we  have  determined 
"he  scaling  conditions  (35).  which  predict  the  final  interval  (in  A  or  in  ~ 1  of  the  scaling. 
Thus,  the  initial  assumption  is  self-consistently  validated.  Taking  into  account  divergence 
at  small  r  in  the  mean-fie'd  equations  (see  the  discussion  after  Eq.  (3)).  the  existence  of  the 
scaling  means  that  in  the  region  (35)  the  large-scale  fluctuations  dominate  in  determining 
the  polarizability  and  cause  failure  of  the  mean-field  approach.  This  situation  is  similar  to 
the  phase  transition  with  X  as  the  transition  parameter.  A'  =  0  being  the  transition  point 
and  the  polarizability  playing  its  usual  role. 

We  now  proceed  to  find  the  relation  between  the  absorption  of  a  fractal.  Im\(A  ), 
and  density  of  its  eigenstates.  u{  A”),  which  is  introduced  in  the  conventional  way  as 

"(A)  =  l£^RY-Wn))  .  (37) 

n 

The  density  v  is  normalized  to  one  monomer,  and  the  mean  value  of  A’  weighted  by  u(X ) 
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is  equal  to  zero  since  the  interaction  tensor  (4)  is  traceless: 

X  50 

J  u\  X )tLY  =  3  .  J  X u{  X )dA'  =  TrU’  =  0  .  (3S) 

—  50  —  OO 

The  value  of  the  first  integral  in  (38)  is  due  to  the  vector  character  of  the  interaction 
between  monomers. 

Isolating  in  Eq.  (28)  the  terms  with  i  =  j  using  the  orthogonality  property  of  the 
matrix  U.  from  Eq.  (8)  and  comparing  the  obtained  results  with  Eq.  (37).  we  obtain 

Im\(  A’)  =  —  t't  X )  +  — —  Y*  !  (n\lr\ia){n\U\ja)  )  .  (39) 

3  3.V  \  :  1  1 

r^j.n 

It  follows  from  (38)  that  the  first  term  in  this  expression  (39)  satisfies  both  sum  rules  for 
Im\,  ( 14)  and  (23).  Correspondingly,  the  zeroth  and  first  moments  in  X  of  the  second  term 
in  (39)  are  exactly  equal  to  zero.  Since  this  term  should  also  have  a  scaling  dependence 
on  A'  (provided  (35)  is  satisfied),  it  then  equals  zero.  Hence,  in  the  collective  region 

Imx(A')  =  ^u(X)  oc  |A’|do-1  .  (40) 

i.e.  both  the  absorption  and  eigenstate  density  scale  with  the  same  index  d0  (the  optical 
spectral  dimension). 

We  emphasize  that  the  relation  (40)  between  the  fractal  optical  absorption  and 
the  eigenstate  density  of  its  dipolar  excitations  has  also  the  character  of  the  theorem  based 
upon  exact  sum  rules  and  the  scaling  assumption.  The  physical  meaning  of  Eq.  (40)  is 
clear.  In  a  strongly-disordered  system  such  as  a  fractal,  all  collective  eigenstates  with 
equal  statistical  weights  contribute  to  the  optical  absorption.  The  relation  (40),  however, 
should  only  hold  for  nontrivial  fractals  in  the  spectral  region  (35).  Note  that  at  the  wings 
of  the  absorption  contour,  which  are  described  in  the  binary  approximation  (Sec.  4),  it  is 
violated  In  ordered  systems,  selection  rules  are  valid,  and  only  few  states  of  corresponding 
symmetry  contribute  to  the  absorption. 

The  scaling  behavior  of  the  density  of  fractal  eigenstates  as  a  function  of  frequency 
was  introduced  by  Alexander  and  Orbach  [16]  for  the  problem  of  mechanical  vibrations 
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of  fractals  ( "fractons"  i.  where  the  index  d .  correspondina  to  <i0  in  (40).  was  called  the 
spectral  (fracton)  dimension.  In  the  general  case,  the  value  of  the  spectral  dimension  is 
determined  by  the  internal  fractal  geometry,  and.  for  a  vector  problem,  also  by  the  manner 
in  which  the  fractal  is  embedded  into  the  three-dimensional  space.  In  principle,  the  spectral 
dimension  may  also  depend  on  the  interaction  between  monomers.  Thus,  the  value  of  the 
optical  spectral  dimension  d  should  not  necessarily  coincide  with  the  vibrational  spectral 
dimension  in  Ref.  (161. 

In  Refs.  [15-17]  different  ways  are  provided  for  finding  an  expression  for  the  dis¬ 
persion  relation  of  fractal  vibrations,  which  in  our  notation  has  the  form 


,-D/J 
A  i  xlv 


(41) 


which  obviously  differs  from  (34).  In  particular,  for  a  trivial  fractal  ( D  =  3).  the  relation 
(41)  predicts  a  gapless  dispersion  law.  which  is  linear  for  d=  3.  An  inspection  of  the 
treatments  of  Refs.  [15-17]  shows  Eq.  (41)  to  be  valid  only  for  excitations  of  the  Goldstone- 
tvpe  (acoustic  phonons,  magnons  and  so  on),  whose  corresponding  creation-annihilation 
operators  in  the  long- wavelength  limit  turn  to  generators  of  symmetry  groups  of  the  given 
system  (translations,  rotations,  etc.).  The  dipolar  excitations  are  not  of  the  Goldstone- 
type.  but  rather  tend  to  surface  piasmons.  which  do  have  a  gap  in  the  spectrum.  In  such 
a  way.  for  the  plasmon-Iike  excitations  considered  in  present  paper  the  dispersion  relation 
(34)  is  valid,  and  not  (41).  which  holds  for  phonon-like  excitations. 

Finally,  we  briefly  consider  Re\(A”).  Since  the  kernel  of  the  integral  (Gilbert) 
transform  (13)  by  itself  does  not  provide  convergence,  both  the  collective  region  (A  <C  1) 
and  wings  of  the  absorption  contour  (A*  >1)  give  rise  to  Re\(A').  Therefore.  Rey(A  )  is 
not  described  by  a  scaling  dependence  even  for  A*  in  the  region  (35).  This  conclusion  can 
also  be  reached  by  starting  directly  from  ( 12). 


4.  Wings  of  the  absorption  contour  and  model  of  diluted  fractals 

In  the  previous  section  the  collective  region  was  considered  in  which  the  eigenstates 
of  a  fractal  are  determined  by  interactions  on  large  scales  and  spread  over  many  monomers. 
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The  observable  quantities  sroverned  by  such  excitations  are  insensitive  to  small-scale  details 
of  fractal  structure  and  the  form  of  the  interaction  of  monomers  at  small  distances.  The 
behavior  of  such  quantities  should  be  universal  and  scaling.  Below  we  shall  describe  the 
wings  of  the  fractal  spectrum  i  the  region  j_V |  >  1  i.  In  this  case  the  external  field  interacts 
with  those  fractal  excitations  which,  in  accord  with  (33).  are  localized  on  a  few  monomers. 
Therefore,  small-scale  details  are  essential. 

It  follows  from  the  above  arguments  that  it  is  in  principle  impossible  to  formulate 
a  universal  description  of  the  wings  of  a  fractal  spectrum.  It  thus  seems  appropriate  to 
employ  a  model  of  a  fractal  with  extremely  simple  structure  on  a  small  scale.  For  such 
a  model  we  introduce  dilute  fractals  which  could  be  obtained  by  roarsing  a  structure  of 
original  fractals.  Instead  of  the  renorm-transformation  described  in  Sec.  3.  we  shall  use 
a  simple  and  numerically  efficient  algorithm  for  the  structure  eoarsing,  called  dilution, 
which  consists  of  the  following.  Each  monomer  of  an  original  fractal  is  randomly  left  at  its 
position  with  some  small  probability  3  <C  1.  or  is  removed  with  the  probability  1  —  3.  Then 
the  fractal  as  a  whole  is  reduced  in  size  (l/f?)1/,£)  times.  After  that,  we  obtain  a  fractal 
which  is  characterized  by  unchanged  critical  exponents  governing  the  collective  region,  with 
the  pair  correlation  function  (1)  (f?o  =  1)  being  valid  up  to  small  distances  on  order  of 
/•i)  =  i  3)l/°  C  1-  The  dilution,  indeed,  does  not  change  the  polarizability  of  a  monomer, 
i.e..  it  conserves  the  variable  Z.  and  the  total  size  of  the  fractal  is  diminished  by  a  factor 
of  r o-  The  dilution  ( random  decimation,  in  other  words)  eliminates  the  local  structure  of  a 
fractal,  such  as  connectivity,  separation  between  nearest  monomers,  coordination  numbers, 
etc.,  without  affecting  the  collective  region  (Sec.  3). 

The  diluted  fractals  do  not  only  provide  a  convenient  theoretical  model,  needed,  in 
particular,  to  check  the  numerical  simulation  in  the  spectral  wings  (see  below).  They  also 
can  be  considered  as  an  idealization  for  a  rather  wide  class  of  fractal  clusters  in  nature, 
which  are  formed  by  random  binding  of  impurities  to  an  original  I  nondiluted)  fractal,  with 
a  small  fraction  of  occupied  binding  sites.  As  examples,  we  may  mention  dye  molecules  or 
metallic  microparticles  in  complexes  with  macromolecules,  or  in  pores  of  glasses,  and  so 
on.  Hence,  the  description  of  spectral  wings  in  the  model  of  diluted  fractals  is  interesting 
by  itself. 
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To  describe  the  interaction  at  small  distances,  which  is  important  for  iA  j  >  1.  we 
shall  use  the  binary  approximation  developed  earlier  iS-llj.  For  each  monomer  we  shall 
exactly  take  into  account  its  interaction  with  its  nearest  monomer  neighbor,  neglecting  the 
effect  of  other  monomers.  The  polarizability  of  a  monomer  in  such  a  pair  is  given  by  the 
tensor  \^j  [8], 


\2  =  (Z  t  6aii  -r  j  (Z  —  2r~i)~l  —  (  Z  4-  r~ 


-i\~l 


-^-■1 


where  rn  is  the  radius  vector  of  the  nearest  neighbor,  n  =  r„ ,rn ,  and  the  superscript  (p) 
stands  for  "pair".  After  averaging  over  r  .  we  obtain  the  mean  polarizability  of  a  monomer 
in  the  fractal  as 

?C 

\=l-J[2{Z  +  r?)-l-(Z-2r?)-1}p{rn)<lrn  .  (43) 

o 

where  pfrn)  is  the  density  of  the  distribution  in  distances  between  nearest  neighbors.  For 
A’  I  »  6  it  follows  from  (43)  that  (cf.  (28)) 


Im\  =  [2p(  rn  )0(  A")  -r  21/<3p(21/3cn  J0(— A" 


!rn  =  A' 


.-1-1/3 


where  0f  A”)  is  the  Heaviside  unit  step  function. 

At  |  A'j  A>  1  the  main  contribution  to  (44)  occurs  at  small  distances  rn  .  Here  the 
following  asymptotic  form,  derived  from  (1).  is  valid. 


p  (  f n  )  ~  Dr®  1 


which  is  universal  for  diluted  fractals.  Substitution  of  (45)  into  (44)  yields  the  following 
expression  for  the  far- wing  absorption  contour: 

Im\  =  —  jA’pl  +  D/3)  r0(A')  +  2O/3"10(— A')|  .  (46) 

9  i 

As  one  can  see.  this  is  a  scaling  form  but.  in  contrast  to  the  collective  region,  the  exponent 
is  greater  than  unity,  and  the  spectrum  is  asymmetric  with  respect  to  the  sign  of  A  (cf. 

(  40)  and  the  discussion  of  it).  One  can  show  by  adding  a  third  monomer  to  a  pair  of 
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nearest  neighbors  that  the  form  of  the  absorption  wing  is  given  correctly  by  the  binary 

approximation  (46)  in  the  leading  term  in  ;.Yj  (the  relative  magnitude  of  corrections  is  on 

,  .  —  2  . 

the  order  of  LY j  "). 

For  the  case  rn  ~  1.  the  distribution  p{rn  j  cannot  be  found  in  a  universal  form 
since  it  is  determined  not  only  by  ( 1 )  but  also  by  higher  correlation  functions.  To  derive 
a  rough  estimate,  let  us  neglect  higher  correlations.  Then  the  number  of  monomers  in  a 
fixed  vicinity  of  a  given  monomer  obeys  the  Poisson  distribution  and.  consequently, 

,nr„  )  =  Dr%~]  exp(-r^)  .  (47) 

Substitution  of  (47)  into  (43)  results  in  the  expression 

Z)  .  1 48) 

where  the  function  Sd  of  the  complex  variable  Z  is  defined  by  means  of  the  integral 
representation 

Sd<Z)  =  -j  J  !/-|l+D,1,exp(-!/-DC)nh-  (49) 

U 

and  its  analytic  continuation.  The  function  Sd{%)  is  analytic  in  the  complex  plane  with 
the  branch  cut  from  0  to  D  along  the  real  axis;  Z  =  0  is  a  second-order  branch  point.  The 
polarizability  (48)  satisfies  all  the  exact  relations  obtained  above. 

The  form  of  the  expression  i48)  coincides  with  that  previously  derived  [S.9],  al¬ 
though  the  function  5q(Z)  is  different.  The  difference  is  due  to  an  incorrect  averaging  in 
[8.9]  over  configurations  of  a  monomer  pair  in  a  fractal  (see  the  text  between  Eqs.  (14) 
and  (15)  in  Ref.  [9]).  \Ve  emphasize  that  all  the  results  of  Refs.  [S- 11]  remain  valid,  except 
that  now  the  form  (49)  of  the  So-function  should  be  used. 

It  can  easily  be  shown  that  in  the  binary  approximation  the  eigenstate  density  is 
symmetric  in  A'  and  has  the  form 

C-Y)  =  \J  [2A(|.V|  -  r-J)  +«(|.V|-2r^)]p(r„Klr„  .  (50) 

The  relation  (40)  between  Im\(A’)  and  i/(.Y)  in  the  binary  region  does  not  hold  unless  the 
polarizability  (43)  is  symmetrized  in  X.  Hence,  the  contribution  from  the  binary  region 
does  not  interfere  with  the  normalization  of  the  relation  (40)  in  the  collective  region. 


,if, 
3  r 


25o(-Z)  -  -5d{  ^ 


21 


5.  Results  of  numerical  simulation 

Three  tvpes  of  fractals  were  studied  as  obtained  by  the  Monte-Carlo  method  with 
the  use  of  a  program  generator  of  random  numbers: 

(i)  Random  walks  ( D  =  2)  were  built  on  a  cubic  lattice  with  monomers  placed 
at  visited  sites  and  subjected  to  1000-fold  ( 3  =  0.001)  dilution  (random  decimation  -  see 
Sec.  4).  The  fractals  obtained  after  the  dilution  contained  a  mean  number  of  monomers 
of  (.V)  =  30  (for  the  control.  (N)  =  15). 

(ii)  Self-avoiding  random  walks  1 1  Si  were  generated  with  fixed  length  and  random 
<tep  direction  usins  a  trial-and-error  method.  Monomers  were  placed  at  nodes  and  taken 
ro  be  hard  spheres  of  diameter  equal  to  the  step.  If  any  step  brought  about  an  intersection 
of  spheres,  then  such  a  step  was  rejected  and  repeated.  A  20-fold  dilution  (  3  =  0.05)  was 
performed.  The  clusters  thus  obtained  contain  a  mean  number  of  monomer  of  (N)  =  25. 
Their  fractal  dimension  D  =  l.S  found  with  the  use  of  Eq.  ( 1)  agreed  reasonably  with  the 
published  value  of  D  —  1.7  [IS]. 

(iii)  The  Witten-Sander  clusters  [19]  were  built  by  the  simulation  of  a  diffusion- 
controlled  aggregation  of  monomers  on  a  cubic  lattice  and  were  20-fold  diluted  (3  =  0.05) 
to  obtain  (.V)  =  25.  The  fractal  dimension  was  found  as  D  =  2.5.  in  close  agreement  with 
the  known  value  of  D  =  2.51  [19]. 

As  the  control,  a  gas  of  monomers  [D  =  3)  was  also  used.  The  dilution  for  cases 
i  ii )  and  (iii)  above  was  not  very  high,  nor  was  the  number  of  monomers  in  the  clusters 
studied  very  large  due  to  limitations  of  existing  computing  facilities.  However,  as  will  be 
shown  below,  even  for  such  fractals  the  polarizability  per  a  monomer  \  does  not  depend 
significantly  on  the  number  of  monomers. 

The  polarizability  was  computed  from  (12)  by  diagonalization  of  the  matrix  (4) 
using  the  Jacobi  method  [201.  The  principal  advantage  of  this  approach  is  that  the  relations 
(  13).  (14)  and  (IS)  are  met  exactly.  Its  merits  as  a  computation  method  are  due  to  the 
fact  that  upon  performing  the  diagonalization  once,  one  finds  the  solution  for  any  A  and 
6  by  simple  summation.  For  the  sake  of  control  we  also  used  a  direct  solution  of  the 
original  system  by  the  Gauss  and  square-root  methods.  For  all  cases  the  three  indicated 


methods  agreed  perfectly.  The  calculated  individual  polarizabilities  were  averaged  over  a 
large  ensemble  of  fractals  I  75  clusters  for  '.VI  =  30  and  i.Y)  =  25.  or  150  for  (.V)  =  15). 
The  statistical  error  achieved  was  less  than  one  percent  of  the  maximum  magnitude. 

The  major  portion  of  our  results  will  be  demonstrated  for  the  random- walk  fractals, 
since  in  this  case  the  theoretical  values  of  the  parameters  are  exactly  known:  D  =  2 
and  Ro  —  R\ (6d)-1 '2 .  where  R\  is  the  lattice  period.  Besides,  such  fractals  can  easily 
be  generated  and  diluted,  and  thus  not  limited  by  computer  capabilities.  Therefore,  a 
high  degree  of  dilution  i  3  =  0.001)  can  be  achieved.  Other  types  of  fractals  will  also  be 
considered,  but  only  to  examine  dependences  on  D  for  the  cases  indicated  explicitly. 

The  iinear  polarizability  \  as  a  function  of  X  is  shown  in  Fiff.  1  for  different  values 
of  the  parameter  f.  One  can  see  that  the  absorption  contour  Im\ i  A" )  has  the  form  of  a  wide 
peak,  with  its  width  decreasing  and  height  growing  with  increase  of  Q  =  1/6.  In  accord 
with  the  results  obtained,  the  dependence  on  for  j.Y|  6  is  weak.  The  obvious  increase 
of  statistical  noise  with  diminishing  6  is  in  agreement  with  Eq.  1,13)  (the  optical  theorem). 
The  width  of  the  fractal  absorption  contour  greatly  exceeds  that  of  an  individual  monomer 
f  ~  d).  The  qualitative  behavior  of  the  polarizability  for  fractals  of  the  other  types  studied 
is  close  to  that  shown  in  Fie.  1. 

The  absorption  contours  calculated  for  the  fractals  with  significantly  different  mean 
numbers  of  monomers  are  shown  in  Fig.  2.  It  can  be  seen  that  for  (.V)  >  15  the  depen¬ 
dence  of  \  on  (.V)  has  practically  leveled  off.  This  behavior  supports  the  idea  of  using  a 
characteristic  such  as  the  polarizability  of  a  fractal  per  one  monomer. 

The  fundamental  prediction  of  the  theory  is  the  universal  scaling  behavior  of  the 
absorption  (30)  and  the  eigenstate  density  (40)  in  the  collective  region  (35).  This  prediction 
is  tested  in  Fig.  3.  From  this  figure  it  follows  that  the  absorption  and  density  of  states,  in 
fact,  obey  Eq.  (40)  with  good  accuracy  expected  from  an  estimate  obtained  with  the  aid 
of  Eq.  ( IS). 

The  values  of  the  optical  spectral  dimension  found  from  Fig.  3  and  similar  data 
for  the  other  types  of  fractals  together  with  corresponding  statistical  uncertainties  and 
values  of  D  are  shown  in  Table  1.  It  follows  from  this  table  that  for  fractals  of  a  given 
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tvpe.  the  indices  of  Im\(A"  i  and  ;  i  A’  dirfer  by  less  than  two  standard  deviations,  which 
:s  statistically  insignificant.  For  ail  the  cases,  the  magnitudes  of  u.,  deviate  significantly 
i  on  the  level  of  five  to  twenty  standard  deviations)  from  the  iimiting  values  oi  </0  =  1  and 
and  </„  =  0.  The  fact  that  d0  =  1  is  evident  directly  from  Figs.  1  and  2:  if  d0  =  1  then 
the  central  portion  of  the  absorption  contour  would  be  net  a  peak  but  constant,  whitn 
is  cieariy  not  the  case.  We  note  the  absence  of  a  significant  correlation  between  d0  and 
D.  For  the  trivial  (D  =  3)  fractal  studied,  no  pronounced  scaling  behavior  is  observed,  in 
corresnondence  with  the  theory  (see.  in  particular,  the  collectivity  criterion  i3ol). 

It  is  shown  above  that  for  diluted  fractals,  the  wines  of  the  absorption  contour 
X  ■  li  shouid  correctly  i  in  the  ieadine  order  in  A  )  be  reproduced  by  the  binary 
approximation.  Reeuiar  computation  of  \i  A  i  in  the  framework  oi  this  approximation  was 
carried  out  according  to  the  following  procedure.  For  each  monomer  the  closest  neighbor 
was  found,  the  quantity  \T'(A" )  i  42)  was  calculated,  then  averaging  over  ail  the  monomers 
of  the  given  cluster  and  the  ensemble  of  clusters  was  carried  out.  The  absorption  found 
'his  way  is  compared  in  Fig.  4  with  the  result  computed  by  the  Monte  Carlo  method 
usine  the  accurate  formula  f  12)  |<*i  =  0.33)  One  can  see  that  the  agreement  is.  indeed,  very 
good  for  A’!  >  2.  Even  at  the  center  of  the  contour  the  differences  are  not  large,  which  is 
due  to  Tie  fact  that  <"  is  not  very  large.  In  fact,  the  maximal  radius  of  coherence  for  the 
above  indicated  dimensions,  which  can  be  estimated  from  (33)  by  substituting  | A  )  =  S. 
is  L  v  ~  2.  and  this  does  not  strongly  exceed  Ro  =  1.  Thus,  in  this  case  the  universal 
-waling  behavior  is  not  developed,  even  at  the  center  of  the  contour,  as  a  consequence  of 
the  violation  of  the  condition  .30). 

In  Fig.  4  the  result  is  also  presented  I  dashed  curve  i  for  a  calculation  of  the  estimate 
provided  by  the  analytic  formula  i4S).  Its  agreement  with  the  result  obtained  from  I  12)  in 
the  wings  of  contour  is  reasonably  good,  although  at  the  center  it  is  worse  than  the  regular 
binary  approximation  i  dotted  curve). 

Similar  data  for  c  =  0.03  are  shown  in  Fig.  5.  Comparing  this  with  Fig.  4. 
we  conclude  that  the  agreement  between  the  three  shown  curves  at  the  wings  (  |A  j  >  2) 
remains  good.  However,  the  central  peak  can  no  longer  be  described  by  any  version  of  the 
binary  approximation.  This  peak  is  symmetric,  wheieas  the  binary  approximation  yields 
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an  asymmetric  contour,  and  given  that  it  is  scalable  I  cf.  Fig.  31.  it  is  certainly  collective. 

Lastly,  let  us  describe  one  of  the  most  important  properties  of  fractals,  namely 
rhe  existence  of  high  local  fields  which  can  sigmricantly  exceed  the  external  field.  Such 
fluctuations  bring  about  giant  enhancement  of  the  Raman  and  nonlinear  parametric  scat¬ 
tering  of  light  from  fractals,  which  has  earlier  been  predicted  based  on  the  use  of  the  binary 
approximation  I S- 111. 

From  the  optical  theorem  i  20 1  and  Eqa.  ( 40 )  and  (461.  for  the  enhancement  fac¬ 
tor  of  the  local  field  acting  upon  a  monomer,  one  can  derive  the  estimates,  which  are 
asymptotically  exact  in  the  sense  of  the  dependence  on  i.Yl  and  ts  in  the  collective  region 


1  1  y  |J,+  1 


and  for  the  diluted  fractals  in  the  binary  region  as 

~  1  ,  „  1 1  —  0/3 

G  ~  -A 
d  '  ! 


We  recall  that  the  existencr  if  a  nonvanishing  collective  region  is  the  property  of  nontrivial 
fractals.  From  the  estimates  (511  and  (52).  it  follows  that  for  nontrivial  fractals  in  both  the 


collective  and  binary  regions,  the  fluctuations  of  the  local  field  are  strong:  the  dispersion 
of  this  field  is  proportional  to  a  large  parameter,  the  O-factor  of  an  isolated  me  nomer  (36). 


The  enhancement  factor  G  as  a  function  of  A'  for  the  tractals  of  the  three  types 
studied  is  shown  in  Fig.  6.  In  accord  with  the  above  presented  theory,  the  curves  in  their 
central  part,  corresponding  to  the  collective  region,  are  svmmetric  and  are  described  by 
rhe  scaling  formula  (51).  The  factor  G  increases  with  |A’|.  in  agreement  with  (51)  and 
‘  52 ) .  reaching  values  on  the  order  of  10  ( for  d  =  0.1)  at  the  wings.  We  point  out  that  there 
exists  no  pronounced  correlation  between  G  and  the  fractal  type,  in  accord  with  the  fact 
rhat  nearly  the  same  values  of  d0  are  found  above  for  different  ’ractals  (see  Table  1). 


In  such  a  wav.  the  present  theory  and  numerical  simulation  results  show  that  large 
fluctuations  of  the  local  fields,  and  hence  giant  enhancement  of  the  Raman  and  parametric 
scattering,  are  not  a  result  of  using  the  binary  approximation  in  [S-lll  but  the  general 
property  of  nontrivial  fractal,  whose  monomers  possess  high-Q  optical  resonances. 
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6.  Concluding  discussion 

In  this  paper  we  described  the  optical  ( dipole)  polarizability  of  fractal  clusters.  The 
principal  question  investigated  is  how  the  seif- similarity  of  a  nontrivial  fractal  in  coordinate 
space  is  transferred  to  the  self-similarity  of  its  dipole  excitations  and  optical  polarizability. 
Our  approach  to  the  problem  of  the  fractal  optical  properties  consists  in  using  the  exa't 
relations  (sum  rules  for  the  polarizability  obtained  in  Sec.  2)  and  consequences  of  the 
scaling.  The  scaling  of  the  optical  polarizability  is  first  introduced  as  an  assumption,  then 
the  scale-invariant  theory  is  developed  (Sec.  3).  Finally,  the  results  of  the  theory  determine 
conditions  of  the  scaling  and  contirm  its  validity.  Xamelv.  the  scaling  behavior  is  shown  to 
exist  in  the  center  of  the  absorption  contour  of  the  fractal.  The  wings  of  this  contour  are 
described  in  the  binary  approximation  (Sec.  4).  The  numerical  simulation  results  (Sec.  5) 
contirm  the  analytical  theory.  The  new  index,  optical  spectral  dimension,  governing  the 
scaling  of  optical  absorption  and.  also,  the  dispersion  law  and  the  density  of  eigenstates  of 
the  dipole  excitations,  is  numerically  determined  for  three  types  of  the  fractals  studied. 

Let  us  discuss  the  principal  points  of  the  present  theory  and  main  results  obtained 
in  some  more  detail,  starting  with  the  basic  system  of  equation  (3).  The  interaction 
between  constituents  of  the  fractal  (called  monomers)  is  described  by  the  dipole-dipole 
interaction  i  see  i  3)  i.  In  the  scaling  region,  where  the  polarizabilities  are  determined  by  the 
interaction  on  distances  much  greater  than  the  monomer  size,  the  dipole-dipole  interaction 
is  universal  and  asymptotically  exact. 

Instead  of  frequency  we  have  introduced  the  new  variable  X  (7).  in  terms 
of  which  the  problem  of  the  optical  responses  is  expressed  in  the  universal  way.  The 
analytical  properties  and  sum  rules  for  the  polarizability  as  a  function  of  X  differ  from 
those  with  u;  as  the  variable.  To  develop  the  theory,  we  introduce  the  decomposition  (see 
Eq.  (12))  of  the  polarizability  in  the  terms  of  the  exact  eigenvectors  and  eigenvalues  of 
the  interaction  operator.  Using  the  analytical  behavior  of  the  polarizability  given  by  this 
decomposition,  tensor  properties  of  the  interaction  and  positiveness  of  absorption  (absence 
of  the  population  inversion),  the  exact  sum  rules  1 14). (IS)  and  (23)  are  established. 

At  this  point,  the  assumption  about  existence  of  the  scaling  behavior  (see  Eqs. 
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>29)  and  (301)  is  invoked.  Using  the  above  indicated  sum  rules  and  the  requirement  of 
the  scale  invariance,  the  forms  of  the  fractal  absorption  i  30 )  and  density  of  the  eigenstates 
i  40)  are  established.  It  is  rigorously  shown  that  both  these  quantities  scaie  with  the  same 
index,  which  is  denoted  as  d0  —  1.  where  d„  is  called  the  optical  spectral  dimension.  The 
index  da  is  similar,  though  it  differs  in  magnitude  from  the  fracton  spectral  dimension  d 
introduced  by  Alexander  and  Orbach  [16]  ( cf. :  both  these  indices  determine  the  density  of 
eigenstates  in  the  same  manner,  but  d0  6  [0. 1],  while  normally  d  £  [1.3]). 

To  make  the  theory  closed,  using  the  found  scaling  properties  (32)  of  A",  we  derive 
the  dispersion  law  for  the  fractal  dipole  excitations  (34).  For  this  derivation,  we  have 
formally  invoked  Alexander’s  strong  localization  hypothesis  that  there  exists  a  single  co¬ 
herence  radius  of  the  collective  excitation.  However,  the  validity  of  this  hypothesis  for  the 
present  problem  follows  directly  from  the  assumption  of  scaling:  only  one  quantity  with  the 
dimension  of  length,  which  obeys  the  required  scaling  properties,  can  be  constructed  from 
all  parameters  of  the  problem.  With  the  aid  of  the  dispersion  law  obtained,  the  condition 
of  the  scaling  for  the  fractal  absorption  is  found  (Eqs.  (35)  and  (36)).  It  is  shown  that  this 
condition  can  be  physically  met  for  the  realistic  model  of  the  monomer  as  a  macroscopic 
sphere.  Thus,  the  existence  of  the  scaling  is  established  in  the  self-consistent  manner. 
More  exactly,  the  scaling  assumption  is  shown  to  be  non-self-contradictive. 

At  the  wings  of  the  absorption  contour  (|A|  >  1).  the  fractal  excitations  are 
primarily  two-particle,  and  the  binary  approximation  [8.9]  is  applicable.  These  excitations 
are  sensitive  to  small-scale  details  of  fractal  structure  and.  therefore,  cannot  be  described 
in  a  universal  form.  In  connection  with  this,  a  model  for  a  diluted  fractal  is  introduced 
which  is  characterized  by  an  extremely  simple  structure  at  small  distances,  completely 
defined  by  the  pair  correlation  (  1).  and  thus  by  the  Haussdorff  dimension  D.  Though  in 
the  binary  region  there  is  no  scaling  in  the  general  case  ( cf.  Eq.( 44 ) ),  the  far  wings  of 
absorption  (see  Eq.  (46))  scales  with  the  index  1  4-  D/3,  which  is  completely  different  from 
the  index  do  —  1  in  the  collective  region  (35).  We  should  emphasize  that  the  behavior  in 
the  binary  region  is  not  universal  and  is  characteristic  of  the  diluted  fractals  (for  other 
fractals,  there  may  be  no  scaling  at  all).  The  scaling  behavior  in  the  collective  region  (35) 
is  completely  universal  and  is  the  same  for  the  diluted  fractals  as  for  the  initial  (before 


dilution  l  ones. 


To  complement  the  theory,  we  have  carried  out  a  numerical  solution  of  the  basic 
equations,  including  averaging  over  an  ensemble  of  fractals.  The  polarizabilities  of  fractals 
over  a  wide  spectral  range  have  been  computed  (see  Fig.  1).  The  central  portions  of  the 
absorption  contour  and  the  eigenstate  density  obey  scaling  with  the  accuracy,  which  is 
expected  from  the  calculation  of  dispersion  (with  the  use  of  Eqs.  (18). (T9)  and  Fig.  6), 
thus  confirming  the  theory.  The  values  of  the  index  d0  are  found  for  the  three  types  of  the 
fractals  studied  (see  Table  1).  The  indices  for  a  given  fractal  found  from  the  absorption 
and  eigenstate  density  do  not  deviate  significantly,  in  accord  with  Eq.  (40).  The  values 
found  for  cf0  differ  strongly  (on  the  level  of  many  standard  deviations)  from  the  limiting 
values  d0  =  0.  1.  Hence,  the  optical  spectral  dimension  d0  proves  to  be  a  new  nontrivial 
index. 

Lastly,  fluctuations  of  the  local  fields  (dipole  momenta)  have  been  computed 
'Fig.  6).  These  fluctuations  proves  to  be  large,  much  greater  than  the  mean  values  of 
the  corresponding  quantities.  This  result  is  in  accord  with  our  understanding  of  the  op¬ 
tical  properties  of  fractals  as  determined  by  strong  fluctuations.  The  scaling  law  and 
magnitude  of  the  fluctuations  in  the  collective  region  are  universal  in  the  sense  that  they 
are  defined  by  d0  and  do  not  depend  on  the  small-scale  structure  of  the  fractal.  The  large 
fluctuations,  in  accord  with  [3-11].  should  bring  about  giant  enhancement  of  the  Raman 
and  nonlinear  parametric  scattering  of  light  from  fractal  clusters,  which  will  be  considered 
elsewhere. 
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Table  1 

The  spectral  dimensions  d0  and  estimates  of  their  statistical  uncertainty  (on  the 
level  of  one  standard  deviation!  found  from  the  absorption  contour  Im\i  A”)  and  eigenstate 
density  u\X)  for  different  fractals  as  indicated.  The  values  of  D  are  given  for  convenience. 


Fractal  type: 
<i0  from  absorption  contour: 
d0  from  eigenstate  density: 

D  found: 


Self- avoiding  walks 

0.43  ±  0.03 


0.54  ±  0.08 


l.S 


Random  walks 
0.38  ±  0.03 
0.33  ±  0.05 
2.0 


Witten-Sander's  clusters 
0.49  ±  0.02 
0.51  ±  0.06 
2.5 


31 


Fin.  1.  Fractal  polarizability  i  D  =  2.  X  =  30)  for  different  parameters. 

Fisc.  2.  Fractal  absorption  \  D  =  2)  for  different  values  of  the  mean  number  of  monomers. 

Fie.  3.  Fractal  absorption  as  a  function  of  A*  (for  X  >0)  'ui  a  double  logarithmic  scale:  1 
-  found  from  the  exact  formula  (12):  2  -  expressed  via  the  state  density  according 
to  (40).  The  corresponding  lines  are  obtained  by  linear  regression.  The  values  of 
the  optical  dimension  obtained  see  in  Table  1. 

Fig.  4.  Fractal  absorption  at  6  =  0.33:  1  -  calculated  from  the  exact  formula  (12):  2  -  in 
the  binary  approximation  using  (42)  and  numerical  averaging:  3  -  in  the  binary 
approximation  from  the  estimate  formula  i4S). 

Fig.  5.  Same  as  Fig.  4  but  for  d  =  0.03. 

Fig.  6.  Enhancement  factor  G  (normalized  by  6  )  of  the  local  field  intensity  plotted  against 
A*  for  the  three  types  of  fractals  studied  (D  —  l.S.  2.0.  2.5)  with  d  =  0.1. 
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